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In this paper, we study the probability that a network within a given geometry is fully connected. We consider two generic 

(N. 

scenarios of practical importance to wireless communications, in which one or more nodes are located outside the convex space 
where the remaining nodes reside. Consequently, conventional approaches with the underlying assumption that only line-of-sight 
(LOS) or direct connections between nodes are possible, fail to provide the correct analysis for the connectivity. We present 
an analytical framework that explicitly considers the effects of reflections from the system boundaries on the full connection 
^ ■ probability. A simple two-dimensional geometry is first considered, followed by a more practical three-dimensional system. The 

ly-^ , analysis presented can potentially be extended to networks residing in a number of non-convex geometries. To corroborate our 
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\Q • derivations, we compare our theoretical results with simulated data. Furthermore, we analyze the effects of different system 
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parameters on the connectivity of the network through simulation and analysis. 



I. Introduction 



Autonomous, self-organizing wireless relay networks have gained a considerable amount of research interest over recent years 
(TJ. The potential of such networks to improve throughput and extend geographical coverage without necessarily requiring an 
increase in transmit power have resulted in them being attractive solutions for many modern wireless systems, including WiMAX 
and LTE-Advanced [2|. With similar objectives, this network architecture has also been proposed for use in communications 
at 60 GHz |3], (4) where severe path loss can significantly limit the transmission range of devices. 

One fundamental requirement of multi-hop relay networks, however, is connectivity, i.e., any node should either be directly 
connected to every other node or via a series of intermediate nodes. Over recent years, a considerable amount of attention 
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has been devoted to understanding the connectivity in such architectures (see for e.g., O— O and references therein). In the 
most simplistic case, the connection probability based on the unit disk model, in which nodes can communicate only if they 
are within a specified distance from each other, has been investigated. In 0, for instance, the critical power such that an ad 
hoc network is asymptotically connected with probability one is derived. The transmitting range for connectivity in sparse 
networks is investigated in J9), while the connection probability in the presence of mobile nodes is presented in |10|. 

On the other hand, more practical models have been studied in |8), lfTTl - lfT6l . In ifTTIl . fl3l . the impact of log-normal 
shadowing on the connectivity is addressed. Means of computing the node isolation probability and coverage for dense ad hoc 
networks in the presence of fast fading is presented in |12|. More recently, the impact of boundaries in dense networks have 
been investigated in [15], JT6), where it is shown that in the dense regime, the geometry of a system should not be neglected 
in analyzing the full connection probability. 

Although the connectivity of random networks have been extensively studied, most existing work has considered only convex 
geometries in which the network resides. In other words, there is always an explicit or implicit assumption that connection 
between two nodes is only possible if a direct path exists between the two. In this paper, we study the scenarios where 
such assumptions need not hold. More specifically, we investigate the connectivity in a network where more than one node 
may reside close to gaps outside of the boundaries. Moreover, considering that in practical situations, wireless connectivity in 
confined geometries may involve reflection of signals from boundaries, we present a means of modeling these reflected rays 
using the principles from mathematical billiards ifTTl , While reflections from other obstacles may occur in practice, these are 
less dominant compared to reflections from large, smooth surfaces such as walls. On that account, these minor reflections are 
not explicitly considered in the derivations presented herein. Rather, they are accounted for by choosing an appropriate channel 
fading model. 

Two scenarios are considered. In the first case, the full connection probability of a network where one node (a receiver) is 
located outside a gap (the so-called 'keyhole phenomenon') while all other transmitter nodes are internal to the main system is 
considered. Under the system assumptions, connection from the outside node to any other node may only be possible through 
reflections. In particular, this would be the case if the size of the gap is small and no line-of-sight (LOS) is present between 
the external and internal nodes. We refer to this problem as the escape problem, analogous to the problem in mathematical 
billiards where a particle escaping a given geometry through a gap in the boundary is analyzed lPT8l . In the second scenario, an 
extension of the escape problem is investigated, where both transmitter and receiver nodes are positioned next to gaps outside 
the main system. This problem is referred to as the transport problem, analogous to the work done in fl9l . From a wireless 
communications perspective, such scenarios are particularly relevant for transmission at high frequencies where communication 
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through walls is likely to render the signal undetectable at the receiver. One potential application is the transmission at 60 
GHz. Nevertheless, it should be borne in mind that the analysis presented herein is applicable to any range of transmission 
frequency and conclusions drawn in this paper may be applicable to other systems as well. 

This paper is structured as follows. In Section [II] we present the general setting for the escape problem and derive the full 
connection probability of the network. We also present some analytical results for this scenario. In Section [Hi] we extend our 
analysis to the transport problem. Finally some concluding remarks are presented in Section [IV] 

II. The Escape Problem 

In this section, we analyze the escape problem. We start by describing the general settings of the system and present the 
means of analyzing the effects of reflections from the boundaries using some of the tools from mathematical billiards. The 
full connection probability is derived next and finally some analytical results are presented. 

A. System Geometry 

Consider a simple system as shown in Fig.Q] where a total of N nodes are randomly distributed within the system boundaries 
and one node is located next to a gap outside the system. Let the width and length be w and L units respectively, where 
we assume that L ^> w and the size of the gap is much smaller than the other system dimensions. For this system, we are 
interested in finding the probability that all the N + 1 nodes are connected to each other. This problem can be split into two 
parts, namely, the probability that all nodes enclosed by the boundaries are fully connected and the probability that the external 
node is connected to at least one other node. While some works have investigated the first part of the problem, the second part 
is still an unexplored issue. Connection between the external node and one other node would be possible only if a signal path 
exists that can 'escape' the main system to connect to the outside node. This problem is similar to the work in [18] where the 
authors investigate the probability that a particle can escape a stadium billiard through a gap along the boundary. 

B. Modeling Reflections 

On the assumption of a very high node density in the enclosed space without an external node, an approach similar to that 
presented in [16] can be followed in the derivation of the full connectivity of the network. However, in situations where few 
nodes are located close to the gap inside the system, or alternately, an LOS transmission is not possible between the external 
and one of the internal nodes, it is more appropriate to study the problem using a different approach. More specifically, the 
contribution of signals reflected from the boundaries should be considered in the connectivity analysis. We aim to demonstrate 
the benefits of considering signal reflections under such conditions. To this end, a proper description of the different reflected 
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paths after a given number of reflections is required. In particular, for a given external node location, the region that can be 
connected given a minimum number of reflections needs to be defined. 

Assume that the gap in the lower boundary of the system, as illustrated in Fig. Q] is small and is far from the either vertical 
boundaries. In this case, only reflections from the horizontal boundaries need to be considered. Given the geometry and the 
non-chaotic behavior of the reflected signals, the method of trajectory unfolding can be employed to determine the reflected 
signal paths ifTTl . The trick in this method consists in reflecting the system (rectangle in our case, c.f., Fig. [TJl instead of 
reflecting the trajectory of the signal. Consequently, the signal path can be viewed as a single, straight line that travels through 
multiple images of the system enclosure. An illustration of this concept is given in Fig. [2] where the external node is referred 
to as node 0, with coordinates (xq, yo). We let the y-coordinate of the lower boundary be zero. In this contribution, we assume 
that node is fixed. From a practical point of view, such an assumption is valid in scenarios where an external device is 
connected to some fixed infrastructure, such as wired systems or to a fixed antenna. 

To determine the regions where connection between the external and internal nodes is possible after a given minimum 
number of reflections, the geometry is analyzed using the polar coordinates system. In this case, every region can be described 
by a range of escape angles from the gap and a range of distances. The escape angle, cj>, is defined as the angle between the 
vertical axis and the signal trajectory. We define the maximum escape angle as 9, as shown in Fig. |2] It should be noted that 
although the representation indicates only the regions to the left of node 0, the process can be readily repeated to cover the 
regions on the right. For the purpose of illustration, however, the region to the right of the external node will not be explicitly 
considered in what follows. 

Although any points within the boundaries can be connected by different paths undergoing different number of reflections, 
we are primarily interested in the strongest path in the analysis of the connectivity. In general, considering the attenuation 
of a signal upon reflection with an object, the strongest path is more likely to be that which undergoes the least number of 
reflections. On that account, we define the various regions within the boundaries based on the minimum number of reflections 
necessary for a signal to reach these. For example, we define the region covered only by LOS paths (no reflections) as T>o, 
while the region covered by one signal reflection, excluding the LOS region, is T>i, For a given number of reflections, c, D c 
is determined by the angle and the distance between the external node and one internal node, r, which are defined as 
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For completeness, we also define the respective regions, D c , in the Cartesian coordinate system. Let Xc represents the point 
of impact of the signal on the horizontal boundaries after i — 1 reflections, when the signal escape the gap with the angle 9 
as shown in Fig. [3] This value can be computed as 

x^ = x - (iw + \y \) tan 9 (6) 

It can be verified that the region to the left of node 0, Do, is bounded horizontally by x = — (y — w) tan (9 + Xc and x — xq 
and vertically by y — and y — w. Similarly, region Di is bounded by x = (y — w) tan 6+Xc and x = — (y — w) tan 9 + Xc , 
while T>2 is bounded by x — — (y — w) tan 9 + x^ and x = (y — w) tan 9 + Xc . The definition of the boundaries of other 
regions follows the same approach. 

C. Probability of Full Connectivity 

For the derivation of the full connection probability, a cluster expansion technique similar to that presented in 1T61 is followed. 

1) First Order Cluster Expansion: We define the set of nodes as S = {0, 1, • • • , N}, where node is the external node and 
is fixed while the remaining nodes are randomly placed within the system. We refer to the probability that two nodes i and 
j, separated by a distance d(ri,ri), being connected as Hy. This pair connectedness probability is dependent on the signal- 
to-noise ratio (SNR), channel gains and other system parameters that will be introduced later. In this work, we use the same 
notations as in [16], where a graph g = {A, L) consists of a set A C S and a collection of direct links L C S A : 

The set of graphs with nodes in A is defined by G , while the set of graphs with nodes in A having largest connected 
component of size j is expressed as G^. 

For any pair of nodes, the trivial probability that they are connected is given by 

1 = Hij + (1 - Hij) (7) 
Multiplying over all links in a set A to determine to probability of all possible combinations leads to 

i= n [tfy+c 1 e ^ (8) 
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Taking the first order approximation of the above expression under the assumption of a high node density as in |[T6l leads to 
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where (I) refers to the average of the quantity / over all configurations. For N randomly distributed nodes with locations 
Yi E V C R 2 , i = 1, • • • , N, this quantity is defined as 



W = \mj v / ( ri ' r 2'"' > r iv) dr x • • • drjv 



(12) 



where V = |V|. 

The last expression in (fTTT i indicates that (^w2p—oYijjt p (1 — Hjp)j} I s ma de up to two terms; the first corresponds to 
the case where all nodes within the system boundaries are connected to each other but not to node 0, while the second term 
represents the case where a node j,j > 1, is not connected to some other node k,k > 1, k ^ j, and node is connected to at 
least one other node. Each of these terms are analyzed below. 

Consider the first average term in (fTTT i. 
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where p — represents the density of the system. On the other hand, 
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Depending on the system configuration, either of the average terms in (fTTT i can have a more pronounced effect on the full 
connection probability. Under the assumptions of a high internal node density and a small gap in the boundary, the probability 
that all nodes being connected would be mostly influenced by whether the external node is connected to at least one other 
node. Consequently, the contribution of ( TBi i on the full connectivity of the network would be negligible. We analyze such 
scenarios in this contribution and focus on the term in ( fl~3b . 

2) Pair-Connected Functions and Mass of Connectivity: Here, the mass of connectivity for the pair connectedness function, 
H 01 , will be derived to aid in the formulation of the full connection probability. Such a function can be viewed as the 
complement of the information outage probability, P ou t, which can be expressed as |20| 

Pout = P(log 2 (l + SNRx l^l 2 ) <R ) 

- "(w^tss 1 ) (15) 

where \h\ 2 is the channel gain between two nodes, i?o is the minimum information rate and the signal-to-noise ratio is a quality 
that is proportional to 

SNR oc G T G R r-^a c (16) 

where Gt and Gr are the antenna gains at the transmit and receive node respectively which equals to unity under isotropic 
radiation, r is the dimensionless distance (relative to the signal wavelength) between the connected nodes, a is the dimensionless 
factor representing the amount of signal attenuation upon collision/reflection with the boundaries, c is the number of reflections 
in the signal trajectory and 77 is the environment dependent path loss exponent. Typical values of r\ is 2 under the free space 
path loss assumption and 77 > 2 for dense or crowded environments |21|. In the previous work analyzing wireless connectivity 
in confined geometries 03], ifTBI . a fair assumption was made that a Rayleigh fading channel is present and the random variable 
\h\ 2 is exponentially distributed. However, since we are also considering reflected signals in this paper, a more appropriate 
model for the channel fading would be the Rician fading model ll20l . Such a fading model is particularly suitable when a 



strong signal can be observed at the receiver along with other weaker signals. The stronger signal typically is the one that 
travels the shortest path, and equivalently, undergoes fewer reflections before reaching the destination node. 
The channel power density of a Rician fading channel is given by ll20l 
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where K is the Rice factor, generally taken to be greater than unity (with typical values around 4 or 5 being reasonable for this 
work), w is a channel dependent parameter and Iq is the modified Bessel function of the first kind. The cumulative distribution 
function, on the other hand, is 
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where Qi(a, b) is the Marcum Q-function of the first kind. It follows that the connection probability for a pair of nodes is 

H = I -P ou t 

= Qx (V2K, y/2(K + l)/3rna- c ) (19) 

where j3 is a small dimensionless constant. 

Let us first consider the connectivity mass function involving the Hqi term. With reference to (fT3l l and assuming a path loss 
exponent of rj = 2 (free space path loss), 

/ Hoidr = V / / '" aX rQi (V2K, y / 2(K + l)/3r 2 a - c ) drdcf) (20) 

To move forward in the analysis and perform the integration, we use the approximation of the Marcum Q— function presented 
in 11221 . i.e., Qi(a, b) is approximated to exp (— e v b^) where v and /i are parameters dependent on the Rice K— factor. It 
follows that 
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where k c = \J2(K + \)(3a~ c , 7(., .) is the lower incomplete gamma function and C is the maximum number of reflections 
considered. Typically, more than five collisions/reflections with obstacles would result in the signal being too heavily attenuated 
to be recovered, especially for high frequency transmission. It should be noted that in d2"Tl >. only the region to the left of node 
in Fig. [2] is considered for conciseness. Extending the results for the general case is straightforward. To enable the integration 
of the lower incomplete gamma function with respect to the parameter (f>, we consider its series expansion. Referring to © 



and under the practical assumption that 6 is small (a 'keyhole' gap), the expansion about = can be taken for the term 
involving rmax- However, considering the definition of r min in (O, it is realized that an expansion about = may not 
be appropriate given that the cosecant function is not defined at the origin. Instead, an expansion about | is taken. The 
corresponding expansion of the lower incomplete gamma functions are as expressed in d22l to d25T l given below. 
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To show the suitability of the expansions for the above integration, we plot the actual and approximated values of the 
corresponding lower incomplete gamma functions in Fig. |4] and [5] It can be observed that over the range of integration 
(c.f., (1211). the approximated values are reasonably close to the actual ones. Similar observations can be made for larger values 
of 9 than those considered in the plot. 

Substituting (|22]i and (J23j into (EB leads to 
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Recall that the first average term in ( fTTb corresponds to the outage probability with respect to node 0, where all internal nodes 
are connected to each other but no connection exist to the external node. We compare the theoretical results in (|26*] | with 



simulation results for this function in Fig. [6] For this plot, values of f3 — 10~ 3 , yo = —2, w = 20, L = 100 and a maximum 
C = 6 reflections were chosen. The simulated results were obtained by finding the number of occurrences over which node 
could not connect to any other nodes while all nodes inside the system were fully connected over 10000 different realizations 
of the channels and node locations. The channel gain between every pair of nodes is modeled following the Rician power 
distribution. The closeness of the theoretical and simulated results validates the derivation of the connectivity mass above. 

At this point, we are interested in investigating the effects of each of the system parameters on this outage probability. We 
first analyze the effects of yo on this value, assuming a = 0.6, /3 = 10~ 3 and K = 4. Results are shown in Fig. [7] As expected, 
a smaller value of \yo\ leads to a decrease in H ° ldri . As node moves closer to the interior, the value of 6 increases, 
resulting in a larger area being covered with successive number of reflections from the boundaries. This in turn increases the 
probability of the external node being connected to one of the internal nodes. Analysis of the plot also indicates that as the 
width of the system is increased (at a fixed density of p = 0.1), the probability that node does not connect to any other nodes 
decreases. Such an observation is attributed to a larger LOS region for systems with larger values of w. Similarly, the areas 
covered by successive reflections are increased as well. As such, there is a higher chance that the external node can connect to 
another node inside the boundaries as the signal travels a shorter path and undergoes a lower number of reflections. On that 
account, the signal attenuation is lower. 

The effects of individual reflections is investigated next. To illustrate the benefits of considering reflections in the connectivity 
analysis, the effect of increasing number of reflections is illustrated in Fig. [8] It can be seen from the figure that the probability of 
node not connecting to an internal node decreases as more reflections are considered, especially for large values of a. However, 
it can also be observed from the plot that the influence of three or more reflections is negligible on this outage probability. 
This observation is due to the value of j3 chosen for the simulations. Choosing much smaller values for this parameter (for 
e.g., 10~ 5 or less) may lead to a higher connection probability as the number of reflections increases. Nevertheless assuming 
a high frequency communication system, a larger value of (3 would arise. Signal quality degradation as a function of distance 
covered is more pronounced. Thus, even though the signal is not attenuated on hitting an obstacle (a = 1), the larger path of 
the reflected ray often imply that the received signal at nodes far from the gap is often unrecoverable. 

Another system parameter expected to affect the connectivity of this network is the gap length. The effect of this parameter 
is investigated in Fig. [9] The smaller the gap size (equivalent to a smaller value of 8), the higher is the probability of node 
not connecting to any other internal node. Recognizing the logarithmic scale of the vertical axis in the plot, it can be deduced 
this outage probability would decay exponentially with increasing size of the gap. 

For the sake of completeness, we present the method of evaluating the second average term in ( TTTb . As previously stated, 
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this average term would be negligible under the system configuration considered herein. This term refers to the probability 
that one internal node is not connected to the remaining nodes and involves the evaluation of the integral J Huydri, where 
Him is the pair-connectedness function between nodes 1 and N, both of which are internal to the system. Since the space 
enclosed by the boundaries is convex in our model, a LOS path or direct connection is possible between any pair of nodeil 
On the account that reflected signals would typically be weaker than the direct or LOS signal, reflections from the boundaries 
are not explicitly considered in this part of the analysis. However, a Rician fading model is still assumed. In that case, 
H 1N = Qi (V2K,ryj2(K + which leads to 



(A (( Xl - x N f + ( Vl - y N ) 2 ) *) d Xl d Vl (27) 
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where A = v(2(K + l)/3)^ 2 . It is realized that evaluating the above integrals can be very challenging, even while using the 
approximation of the Marcum Q— function as performed earlier. While the connectivity mass under the assumption of Rayleigh 
fading for a square has been derived in ||T5l . it is not possible to directly extend the work to the Rician fading case. Given the 
difficulty in evaluating the integral of the above exponential term, another approximation to the Marcum-Q function is made, 
namely, Q(a, b) rj exp (— e V2 b 2 \. Compared to the previous approximation, the exponent on the variable b has been changed 
from /i to 2. Similar to the method used in [22J, the derivation of the parameter V2 for a given value of K can be obtained 
by minimizing the sum of square error (SSE) between the actual and approximate function. With this approximation, 
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where f(x, I) — (erf ((I- x)V\ \ + erf (xVx\\. The first integral term in ( TBI ) then becomes 
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To evaluate the above integral, the series expansion of J(xn, L) and f(yN, w) about L/2 and w/2 respectively can be taken, 
which corresponds to the points at which the functions are maximum. 

f(x N , L) » 2erf (^^) - (*N - ^) n + (x N - ^\ 

f(y N , w) » 2erf (7a|) - (y N - |) ' r 2 + (y N - |) ' (30) 

'Throughout this paper, we have assumed that any node does not shadow the signals from other nodes. 
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where t\ — -^L\^e A « and t 2 = ^wA^e X ™ . Neglecting the terms with the product of variables xn and yjy in the 
expansion of f(xN,L)f(yN,w) leads to 
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where ai = p^^-erf ^a/A-j^, a 2 — pr 2 ^erf (^/X^fj and -& — tan 1 (^j-^=j. The last expression in (|3TT l is obtained by 
changing the integration to the polar coordinate system. Using 
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-03 + (| - 0) 



A 2 tf- 1 



(34) 



We now look at the last term in (TBI ), p J J Hqnc ~ p f HlNdri dxNdyN- This integral involves two terms, namely Hqn and 
Hin, which corresponds to the probability of the external node connecting to an internal node and the probability of any pair 
of internal nodes being connected. However, as we have identified above, reflections from the system boundaries are considered 
only for the H n term. Consequently, the integral comprises of a sum of integrals over different D c , for c = 0, 1, • • ■ , C as 



p I Hon ( 1 - y I H 1N dr 1 



N-l c . 

dr N =pJ2 e-~ x < f e- p S HlNdri dv N 



c=0 



(35) 



'D ,, 



where f is the effective distance between node and a corresponding internal node and A c = 2{K + l)(3a c v. Based on our 
previous observation, we consider a maximum of 2 reflections and assume that node lies in the middle of the gap. Assuming 
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that T> c is bounded in the a;— direction by x = and x = (c.f. Section II), the required integral becomes 



p I H 0N ( 1 - — / H 1N dri ) dr N = 



N-l 



2e 



-pf erf(VA-|)erf(\/Af ) 



j u b° e -\){{xN-x ) 2 + (y N -y ) 2 ) e ai{y N -^ Y ^2 (x N ~ ±) 2 fa ^ 
h 

_|_ J u b e -^l((xN-x„) 2 + (2w-y N ~y ) 2 ) e a 1 [y N -f) 2 e a2(x N --^Yfa^ 
i f u b 2) e -*2((x N -xo) 2 + (2w+y N ~y ) 2 ) e a 1 (y N -f) 2 e a2(x N -±) 2 ( i x 



The first term is the square bracket in $36l can be expressed a 



dy N (36) 



o Ji 



-*o((xN-x ) 2 + {y N -y ) 2 ) e <j 1 (y N ~§y e <j 2 (x N -±) 2 fa^^y 



N 



- e -1*-1l 



Jl 



e (o- 2 -Ao)(2;jv-Px) 2 + (o-i-Ao)(yiv-pJ, 0) ) 2 ^ a;Ar ^y 



N 



a/ ((Ti - A )(cr 2 - Aq) J-y/<T 1 -\ p { y ) "'V CT 2-A (;< 0) -p x ) 



(37) 



where 



Px — (a 2 -\ ) ' 9:1 — 20 ( CT2 _a ) ' P» — (cti-Ao) * J — 10 (<ti-a ) ' smmal expressions can De aenvea 



(xQ-L/2) 2 (0) _ ^tu-Aoyo 



(0) 



for the other terms in d36*b . 



Z). Extension to the 3-D Case 

In the previous subsection, a simple 2-D system was considered. Nevertheless, extension of the analysis to more practical 
3-D systems is straightforward, as will be shown here. Consider a system model as shown in Fig. [I0]and let the gap in the 
lower boundary take any shape. Similar to the 2-D case, we assume that the length and depth of the system, L, is much greater 
than the height w for the same reasons as previously mentioned. 

In this subsection, we focus primarily on the probability of the external node being connected at least one internal node (c.f., 
the first average term in (fTTTi). To define the regions within the system boundaries covered by a given number of reflections 
from the boundaries, we make use of the spherical coordinate system such that any point within the geometry can be defined 
by three variables, namely, r$, <f> and ip. The first variable, r^, represents the distance from the external node to an arbitrary 
point inside the geometry; <f>, in this case, corresponds to the latitudinal displacement with respect to the ir— axis while the 
longitudinal displacement with respect to the z— axis is defined by the angle ip. The latter variable is dependent on the shape of 
the gap and also defines the maximum value that <j> can take, which is 6(ip). Considering a fixed external node at (xo,yo, zq), 

2 Since we are not interested in the average term in 1141 in this contribution and due to space restriction, we do not present the full simplification of this 
expression herein. 
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any point (in Euclidean space) inside the boundaries can be expressed as a function of this 3-tuple as 



x(r,<j>,(p) = x + cos ((f) sin (\4>\) 
V{r,(j>,ip) = y + r> sin (<,£>) sin (|0|) 
z(r,</>,<p) = z cos (|0|) 



(38) 



Analogous to the 2-D case, the region D c for c reflections is defined by the following range of distances: 



J«0 
min 



-(c) 



[go] 

COS (ft ' 

2(e«;+|2 |)sm(6>(y)) 
sin(0(»+<« i 

{c+l)w + \z \ 
cos(0) 



c = 



c> 



(39) 



while 



mm 



o. 



c = 



(40) 



tan -l ((e-l)«)+|zo|)tan(8(y)) Q 
(c+l)tu+|zol 



Assuming a Rician fading channel as in the 2-D scenario, the connectivity mass (c.f., $2l[ ) can be evaluated by the following 
integral: 

C r 27T r 8( v ) r r£l 



H 01 dr ~ y2 / / e^q>(~e u Kc^r^r 2 sin (/) dr dcj) dip 

„_n JO J(i (c) Jr {c) 



c=0 

E / / — 



- ( ~, A c (r max ) M ^ - 7 f ~, A c (r min ) M 



sin </> d0 



(41) 



where A c = e"^. Using similar approximations as in ( |23T l and assuming a circular gap on the lower face of the system, the 
integral evaluates to 



Hnidr w 2vr-^ 



A n * 



- cos (%)) ( 7 (f , Xo(w + \z \n) - 7 (f , Aoko|) M ) 
+ (-2 - (-2 + %) 2 ) cos(%)) + 2%) sin(%))) § a| 
x ((tu + \z \) 2 e- Xa ^ Z0 ^ - |z | 2 e- Ao|ao1 ") 

(7 (}, A c ((c + l)w + |*|)") - 7 (j, 2»\ c (cw + \MTi$8$) ' 
x (cos(^t) - «*(%))) 

+f a| ((c + l)w + \z \) 2 e -M(c+i)™+\Mrc 3d + cot(30&)/2) A 3d D 3d + A 3d B 3d E 3d 



(42) 
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where 



C 3d = - {6{v? - 2) cos (%)) + ((^t) 2 - 2) cos (^1) + 2%) sin (%)) - 2^ n sin (> 
Aw = -' 



min 



— 



^ cos (*(„)) + - cos + sin (*(„)) - sin ) 

i (9{ v f - 8) cos + ~ (-8 + (%) - 2<4t)) cos (^ n ) + %) sin (%)) 

+ (%) - 2<f><±) S in(^ n ) (43) 

Using a similar procedure as in the preceding subsection, we analyze the suitability of the derived expression by comparing it 
with results from Monte Carlo simulations. Results are shown in Fig.QT] In this plot, a circular gap of radius 0.1 is assumed, 
zq is set to —2, (3 = 10~ 3 and K = 4. Similar to the 2-D scenario, the plot indicates that the approximations used in the 
derivations of the connectivity mass is appropriate for the analysis of the system. We next analyze the effects of different 
parameters on the investigated outage probability. 

Following the same approach as for the 2-D scenario, we analyze the effects of the reflected rays on the outage probability 
of the system. This effect is illustrated in Fig. [12] where the same parameter values as above are used. As previously observed, 
the contributions of signals undergoing more than 2 reflections on the connectivity is negligible given the system parameters 
(c.f., Fig. EQ. 

On the other hand, the effect of the radius of the gap is illustrated in Fig. [13] while the relation between the outage probability 
and the distance of the external node from the gap is shown in Fig. [14] Similar to the 2-D system, an increase in the value 
of w (corresponding to the height of the system) or increasing size of the gap results in a lower outage probability. The 
reasons for such behavior is as previously stated, namely, a larger LOS region. Nonetheless, the rate of decrease of the outage 
probability in the 3-D case differs from the simple 2-D case. Such an observation can be credited to the way in which the 
area (or volume) of the space covered by signals undergoing reflections with the boundaries behave in the two systems. To 
illustrate this point, let us consider the ratio of areas (or volumes) covered by the LOS signals and signals that are reflected 
once from the boundaries. For the 2-D case, we refer to the area covered by the LOS signals by ao while the area covered 
after one reflection is referred to as a\. Similarly, we use vq and vi to represent to respective volume of space in the 3-D 
system. It can easily be shown that 

2± = (2w+ \yp\f ~ (w+\y a \f 
a o (w + \yo\f -\yo\ 2 
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while 

vi = {2w+\z \Y-(w + \z a \Y 
vo (w + \z \) 3 - \z Q \ 3 

Clearly, the rates at which the connection probability changes while considering reflections are different in the two investigated 
scenarios. 

It should be noted that in the system model under investigation, there are numerous parameters that affect the full connectivity 
of the network. Due to space restrictions, it is not possible to delve in each possible combination that would increase the 
connection probability. The reader should however bear in mind that further analysis can be readily carried based on the 
expressions presented in this paper. 

III. Transport Problem 

In this section, we study the connectivity in the second scenario, namely when both transmitting and receiving nodes are 
located outside the main system. This study is comparable to the work done in [19| where the transport problem is investigated 
for the stadium billiard. Such a scenario can have a number of practical applications. For instance, it may model the situation 
where nodes are located in different rooms next to keyholes. While it is possible to consider the case where multiple nodes are 
located within the main system and derive the full connection probability as for the escape problem, we restrict our analysis 
to the case where nodes are only located outside the system. In doing so, we avoid the repetition of the results presented in 
the preceding section. The derivation of the full connection probability of all nodes within the boundaries, and those external 
to the system close to the two gaps is a simple extension of the escape problem presented above. 

We investigate two cases for this problem. In the first case, it is assumed that the gaps are on opposite sides of the boundaries, 
while in the second case, the gaps are assumed to be on the same side. In either case, the consideration of reflected rays has 
a considerable impact on the connectivity of the network. 

A. System Geometry 

1) Case 1: Gaps on Opposite Side: Consider the system model shown in Fig. [15] Let the two nodes located at positions 
{ x OiVo) an d { x iiDi) represent a receiver and transmitter node respectively. Assume that the positions of the gaps are fixed 
with the corresponding coordinates on the lower boundary of xi x and xi 2 and x Ul and x U2 for those on the upper boundary. 

Depending on the locations of the gaps and the nodes, a LOS transmission may not be possible, in which case, the only 

way for the nodes to communicate would be through reflections from the system boundaries. Similar to the escape problem, 

we assume that X\ x < Xq < xi 2 and x ri < x\ < x U2 and let the maximum angle of a ray escaping the gap measured from the 

vertical axis be 9 = tan -1 x ", V 1 . 

\yo\ 
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Referring to the figure, it is obvious that a direct connection is not possible if xq — (\yo\ + w)tan9 > x U2 . Under such 
conditions, the only trajectories for the signal would be through an even number of collisions. Assuming the signal undergoes 
a total of c reflections, the range of angles for a signal escaping the lower gap and entering the upper gap is 

= tan U+i)^+toJ l46) 

At = mi n(Man-( (c+ y; |>o| )) ,47, 

The corresponding set of distances traveled would then be 

(c) _ (c + l)w + \y \ 



COS ( 



(48) 



(C ) *o - g "i (49) 

It should be noted that in the above, it was assumed that the gaps are not exactly opposite each other. If this was the case, 
only direct rays between the nodes would be possible since no reflected rays would enter the upper boundaries after an even 
number of reflections. 

2) Case 2: Gaps on Same Side: In the second scenario, it is assumed that the gaps are on the same side of the boundaries. 
Thus, connection between nodes outside each gap is only possible after an odd number of reflections as shown in Fig. [16] 
Let the positions of the gap be Xi t , xi 2 , xi 3 and xi 4 . Let the positions of the transmitter and receiver nodes be (xo,2/o) and 
(xi,y±) respectively. For a given transmitter node location, the maximum escape angle from the gap, 9, is given by: 

= tan- 1 ( ShLZ^A (50) 



l2/o| 

The minimum escape angle which would allow more than one reflection to occur with the boundaries is: 

min =tan- 1 ( p^r ) (51) 



2w + | j/o 

Given such a geometry, the connection between the two nodes would be possible only after an odd number of reflections. The 
angle (f^ c \ measured from the vertical axis, for the transmitter and receiver node to connect after c reflections is given by 

0W = tan- 1 ( 7 ' 1 -, 8 °, , ] (52) 



(c + l)w+ \yo\ + | yi 

where c is odd. The minimum number of reflections such that connection would be possible satisfies the relationship that 
0min < <f>^ < 6. For each such path, the distance between the transmitter and receiver nodes is 

r (c) = [(c+l)™ + |3/i| + |3/o|]sec^ c > (53) 

Alternately, the same approach as in the previous section can be followed to determine the region covered in the receiver 
gap for a given transmitter node position and a given number of reflections. In this case, the minimum and maximum escape 
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angles of a signal from the transmitter gap are given by: 



<4 C L = min^tan- 1 )" ^ *° )) (55) 
V \(c+l)w + \y \J J 



while the range of distance for a given angle <j> is 

( C ) x h - x _ (c + l)w+ Ij/oj 



sin </> cos ( 

^(c) _ : r; 4 - Xp 
sin0 



(56) 



where xu, xi 4 and xo are assumed fixed. 



B. Analysis of Connection Probability 

Let a single node be located close to each gap. Given the geometries considered, it can fairly be assumed that a Rician 
fading channel is present, where the signal undergoing the shortest path and minimum number of reflections is the strongest 
(equivalent to a LOS) signal. Similar to the derivations in the previous sections, the connection probability between a transmitter 
and receiver node given a minimum acceptable data rate of i? is i?oi ~ ex P (e _!y 6 A1 ). The connection probability averaged 
over all configurations between the pair of transmitter and receiver nodes can be expressed as: 

[ [ Hmdv^vo (58) 
Vo Vi J ro J Tl 

where Vo and V\ are the areas of the regions where the transmitter and receiver nodes can reside. It is straight forward to 
extend the above derivations to determine the full connectivity if nodes are present within the system boundaries. In that case, 
each gap can be analysed using the same approach as the escape problem as derived previously. 

1) Case 1: Assuming that the gap locations and sizes are fixed, for a given node position, the integral of H m yields 

H idr 1 = \ / / rH 01 drd<f) 

c is even 
c is even 
c is even 

where c m j n is the minimum number of reflections that would allow connection between the two nodes. To perform the integration 
with respect to <f>, a series expansion can be used again. However, in contrast to the escape problem where it was reasonable 
to assume that 4> would be close to 0, such an assumption may not always be valid in this scenario. At the expense of less 
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x 



sine 



(c+l)w + \y \ 



COS( 



d<p (59) 



concise expressions, the expansion is performed about 



leading to: 



•2 fxo^\\ 
fj, \ sm0 



7( -, A c ((c +l)w + \yo\)Becr(<f>) 



where 



7 ^-,csc^(0 c )(x o - x Ul Y^cj - A t ((f> - c )cot(0 c ) 
+ ^(0-0 c ) 2 + O(^-0 c ) 3 



7 I -, A c ((c + l)w + \yo\f sec^(0 c ) ) + D t tan(0 c )(0 - 
DtEt (a 1 ^2 



B t = csc 2 (0 c ) 2 + cos(20 c ) - ^A c cos 2 (0 c )(iro - x Ul ) M csc M (0 c ) 
A = e-^ac+iV+lyoirsec^^t (( c+ i) w + | yo |)2 S ec 2 (^ c ) 
E, = s<r-|«.,.i -2 + cos(20 c ) +/iA c ((c+ l)w+ |j/ o |) A 'sec^(0 c )sin 2 (0 c ) 



(60) 



Substituting the above approximations in d59l ) results in: 



/ H 01 dv 1 = - V ( 



max Tmin 



c-c n 
c is even 



- ( ^,csc^(0 c )(a; o -a; ul )^A c ^ - 7 A c ((c + l)w + \y \f sec^(0 c ) 



( A t cot(^ c ) + A tan 
1 



/.(c) >>2 
''max J 



^rnax^c 



mm 



1 _L wl,( C ) I 



6 



(A(B t + AA) 



(41 



+ 0(41 



<h {c) . ) 4 

max T mm ; 



1 C 

r C— C m i n 

c is even 



.w -<a (c ? 

max Vmin 



(61) 



The novelty in this approach to analyzing wireless connectivity is the consideration of signal reflections from the boundaries. 
To illustrate this effect, the contribution of each received signal to J Hoidri after undergoing a given number of reflections is 
shown in Fig. [TTjfor different values of w. For these plots, a Rician fading model with K = 5 was assumed, with f3 = 10~ 3 , 
a = 0.85, \yo\ — 2, gap lengths of 0.3, x\ x — 15, x Ux — 14.5, while the corresponding value of 9 = 0.0749. Under these 
simulation parameters, the minimum number of reflections for every value of w considered is 0. Thus, the reference to the 
first received signal in the plot corresponds to the LOS signal, the second received signal is the one that undergoes reflections 
and so on. As expected, the effect of signals undergoing higher number of reflections is less significant. 

We next investigate the effect of the size of the gap on the connectivity mass. This is illustrated in Fig. Q~8] Similar to 
previous observations, it can be seen that the larger the gap, the larger is the value of the integral J Hoidri. Consequently, the 
outage probability decreases with increasing gap length. Such a behavior is intuitive since a larger gap length is synonymous 
of a larger LOS region. 
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In Fig. [19] we present the effect of the distance yg on the connectivity mass of the system. For these sets of simulations, 
we assume that a — 0.85. From the plot, it can be observed that as node moves closer to the boundary of the system, the 
connectivity mass decreases, which implies a larger outage probability. Such an observation is different to what was observed 
for the escape problem. This behavior can however be easily explained by considering the geometry of the system. Given 
the positioning of the gaps (same as in the above simulations), moving node closer to the gap actually leads to smaller 
coverage area for the LOS region which leads to a smaller connectivity mass. On the other hand, moving this node closer to 
the boundary also leads to a wider area being covered by the reflected paths that can enter the gap on the upper boundary. 
This explains the different rates at which the two curves decrease in the figure. 

IV. Conclusion 

In this paper, the connectivity of networks in confined geometries has been investigated wherein signal reflections from 
the boundaries have been explicitly considered. In particular, the full connection probability of a network where at least one 
node is located outside the main geometry has been investigated. This paper has provided an in depth mathematical analysis 
of how the connectivity of such networks is affected by various system parameters, such as the physical dimensions of the 
geometry, the size of the opening on the boundary and the reflectivity properties of the boundaries. Although the connectivity 
was initially studied for a 2-D system, it was shown that the analysis can be easily extended to more practical 3-D systems. 
The expressions derived can be calculated to arbitrary precision using standard libraries. Further work on the simplifications 
of the expressions provided herein are ongoing. It is believed that this contribution will provide the framework for further 
analysis of the connection probability of networks residing in convex as well as non-convex geometries. 
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Fig. 1. System model for the 'escape problem', where one external node needs to connect to at least one internal node. 



20 



References 



[1] M. Haenggi, J. G. Andrews, F. Baccelli, O. Dousse, and M. Franceschetti, "Stochastic Geometry and Random Graphs for the Analysis and Design of 

Wireless Networks" IEEE J. Select. Areas Commun., vol. 27, no. 7, pp. 1029-1046, 2009. 
[2] W. Fu, Z. Tao, J. Zhang, and D. P, Agrawal, "Error Control Strategies for WiMAX Multi-Hop Relay Networks," in Proc. IEEE Global Telecommunications 

Conf. GLOBECOM 2009, 2009, pp. 1-6. 
[3] J. Qiao, L. X. Cai, and X. Shen, "Multi-Hop Concurrent Transmission in Millimeter Wave WPANs with Directional Antenna," in Proc. IEEE Int 

Communications (ICC) Conf, 2010, pp. 1-5. 
[4] J. Qiao, L. X. Cai, X. S. Shen, and J. W. Mark, "Enabling Multi-Hop Concurrent Transmissions in 60 GHz Wireless Personal Area Networks," IEEE 

Trans. Wireless Commun., vol. 10, no. 11, pp. 3824-3833, 2011. 
[5] P. Gupta and P. R. Kumar, "Critical Power for Asymptotic Connectivity," in Proc. 37th IEEE Conf. Decision and Control, vol. 1, 1998, pp. 1106-1110. 
[6] Y.-W. Hong and A. Scaglione, "Critical Power for Connectivity with Cooperative Transmission in Wireless Ad Hoc Sensor Networks," in Proc. IEEE 

Workshop Statistical Signal Processing, 2003, pp. 13-16. 
[7] P. Li, C. Zhang, and Y. Fang, "Asymptotic Connectivity in Wireless Ad Hoc Networks Using Directional Antennas," IEEE/ACM Trans. Networking, 

vol. 17, no. 4, pp. 1106-1117, 2009. 

[8] G. Mao and B. D. O. Anderson, "Towards a Better Understanding of Large-Scale Network Models," IEEE/ACM Trans. Networking, vol. 20, no. 2, pp. 
408-421, 2012. 

[9] P. Santi and D. M. Blough, "The critical transmitting range for connectivity in sparse wireless ad hoc networks," vol. 2, no. 1, pp. 25-39, 2003. 
[10] P. Santi, "The Critical Transmitting Range for Connectivity in Mobile Ad Hoc Networks," vol. 4, no. 3, pp. 310-317, 2005. 

[11] C. Bettstetter and C. Hartmann, "Connectivity of Wireless Multihop Networks in a Shadow Fading Environment," Wirel. Netw., vol. 11, no. 5, pp. 

571-579, Sept. 2005. [Online]. Available: pittp://dx.doi.org/10.1007/sll276-005-3513^x1 
[12] D. Miorandi, "The Impact of Channel Randomness on Coverage and Connectivity of Ad Hoc and Sensor Networks," IEEE Trans. Wireless Commun., 

vol. 7, no. 3, pp. 1062-1072, 2008. 

[13] P. Stuedi, O. Chinellato, and G. Alonso, "Connectivity in the Presence of Shadowing in 802. 1 1 Ad Hoc Networks," in Proc. IEEE Wireless Communications 

and Networking Conf, vol. 4, 2005, pp. 2225-2230. 
[14] D. Miorandi and E. Altman, "Coverage and Connectivity of Ad Hoc Networks in the Presence of Channel Randomness," in Proc. IEEE 24th Annual 

Joint Conf. of the IEEE Computer and Communications Societies INFOCOM 2005, vol. 1, 2005, pp. 491-502. 
[15] J. Coon, C. P. Dettmann, and O. Georgiou, "Full Connectivity: Corners, edges and faces," Jan. 2012. [Online]. Available: http://arxiv.org/abs/1201.3123 
[16] J. P. Coon, C. P. Dettmann, and O. Georgiou, "Connectivity of Confined Dense Networks: Boundary Effects and Scaling Laws," vol. arXiv:1201.4013, 

2012. [Online]. Available: |http://arxiv.org/abs/1201.4013vl | 
[17] N. Chernov and R. Markarian, "Chaotic Billiards," in Mathematical Surveys and Monographs, 127. American Mathematical Society, 2006. 
[18] C. P. Dettmann and O. Georgiou, "Survival Probability for the Stadium Billiard," Physica D: Nonlinear Phenomena, Oct. 2009. [Online]. Available: 

|http://dx.doi.org/10.1016/j.physd.2009.09.019| 
[19] , "Transmission and Reflection in the Stadium Billiard: Time-Dependent Asymmetric Transport." Phys Rev E Stat Nonlin Soft Matter Phys, vol. 83, 

no. 3 Pt 2, p. 036212, Mar 2011. 
[20] J. Proakis, Digital Communications, 4th ed. McGraw-Hill Higher Education, Sept 2000, iSBN-13: 978-0072321111. 

[21] J. P. Coon and C. P. Dettmann, "On the Connectivity of 2-D Random Networks with Anisotropically Radiating Nodes," IEEE Communications Letters 
(to appear). 

21 



[22] M. Bocus, C. Dettmann, and J. Coon, "An Approximation of the First Order Marcum Q-Function with Application to Network Connectivity Analysis," 
ArXiv e-prints, vol. 1210.2067v2, Oct. 2012. [Online]. Available: |http://arxiv.org/abs/1210.2067v2 1 



Image of node 













i 







Fig. 2. Illustration of how the method of unfolding can be used to determine the regions of interest for a given number of reflections. 




Fig. 3. Illustration of the regions covered by different number of signal reflections with the boundaries. 
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K = i,B= l(r 3 , e = 0.0831 



■< 1 .5479 



- Exact Value 

- Approximate Value 
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Fig. 4. Comparing 7 \ \i -^c' r maxj with the approximation in {22} for c = 1. 



K = 4, 0= 10- 3 , 9 = 0.0831 
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Fig. 5. Comparing 7 Mp Acr^j^j with the approximation in {23} for c = 1. 



Gap length = 0.3, m = 20. L = 100, y = -2, p = 0.1 
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Fig. 6. Comparison of theoretical and simulation results for the probability of the external node not connected to any internal node, C = 6. 



23 




Fig. 7. Effect of yo on the probability of node not connecting to any internal node. 




Fig. 8. Effects of the number of reflection on the probability of node not connecting to an internal node, w = 15, yo = —2, L = 100, f3 = 10 3 and 
p = 0.1. 
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p = 0.1, L = 100, w = 20 




■ B - LOS only 

■ O - LOS + 1 reflection 

■ A - LOS + 2 reflections 



1.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Gap Length 



Fig. 9. Effect of gap length on the probability of the node not connecting to an internal node. 




( X 0' ^0' Z o) 



Fig. 10. System model for 3-D case; the gap in the lower boundary can take any shape. 



p = 0.08, L = 100, w = 20, zo = -2, Gap radius= 0.1 




Fig. 11. Comparison of outage probability based on theoretical results and simulation data for the 3-D case, with w = 20, L = 100, /3 = 10" 3 and 
2 = -2. 
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Fig. 12. Analysis of the contribution of reflected signals on the outage probability for the 3-D case. 



p = 0.1, L = 100, w = 20, zo = -2 




Gap Radius 



Fig. 13. Effect of the size of the gap on the outage probability for the 3-D case. 




Fig. 14. Effect of the distance of the external node to the gap, Zo, on the outage probability in the 3-D system. 
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7j l 2 

Fig. 15. System model for the transport problem; case 1, where the gaps are on opposite sides of the boundaries. 




w = 10, 5>o = -1, j8 = 1(T 3 , if = 4 



- LOS coimection only 

- LOS + 2 reflections 

- LOS + 4 reflections 




Fig. 18. Effect of gap length on the connectivity mass for the transport problem. 



GapLength = 0.3, w = 10, 13 = 10" 3 , A" - 4 



- LOS connection only 

- LOS + 2 reflections 




Fig. 19. Effect of yo on the connectivity mass for the transport problem. 
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